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Abstract 

We present spectral calculations of nuclear matter properties including three-body forces. Within 
the in-medium T-matrix approach, implemented with the CD-Bonn and Nijmegen [l|Q] potentials 
plus the three- nucleon Urbana interaction , we compute the energy per particle in symmetric and 
neutron matter. The three-body forces are included via an effective density dependent two-body 
force in the in-medium T-matrix equations. After fine tuning the parameters of the three-body 
force to reproduce the phenomenological saturation point in symmetric nuclear matter, we calculate 
the incompressibility and the energy per particle in neutron matter. We find a soft equation of 
state in symmetric nuclear matter but a relatively large value of the symmetry energy. We study 
the the influence of the three-body forces on the single-particle properties. For symmetric matter 
the spectral function is broadened at all momenta and all densities, while an opposite effect is found 
for the case of neutrons only. Noticeable modification of the spectral functions are realized only 
for densities above the saturation density. The modifications of the self-energy and the effective 
mass are not very large and appear to be strongly suppressed above the Fermi momentum. 

PACS numbers: 21.30.Fe, 21.65.-f, 24.10. Cn, 26.60.Kp 
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I. INTRODUCTION 



Many-body calculations based on realistic bare nucleon-nucleon potentials are able to re- 
produce qualitatively but not quantitatively the saturation properties of symmetric nuclear 
matter. The saturation points calculated with different approaches implemented with vari- 
ous choices of the nucleon-nucleon (NN) potential happen to lie on the so called Coester line 
in the energy per particle-density plane, away from the experimentally allowed values {4]. 
The theoretical predictions give a saturation density sensibly higher than the experimental 
value po ~ 0.16 fm~ 3 (usually in the range 1.5 po — 2 po) an d often overbind the nuclear 
system (up to 25%), failing to get close to the empirical binding energy E « — 16MeV. 
These discrepancies and uncertainties get amplified when calculating the equation of state 
(EOS) of pure neutron matter, which is necessary for the estimates of key quantities such 
as the symmetry energy and in general for the description of neutron-rich matter in neutron 
stars. 

It is believed that these deficiencies originate from neglecting three-body forces (TBF) 
between nucleons in the dense medium. Variational and Brueckner-Hartree-Fock (BHF) 
calculations which take TBF into account achieve a realistic description of cold nuclear 



matter, reproducing: t 



accuracy 



re saturation point of the symmetric EOS with a good degree of 



10 
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12. 



13|- 



The self-consistent real-time Green's functions technique can be used practically to per- 
form calculations of in-medium propagators and properties of a many-body system. In 
nuclear matter this approach considers the in-medium T-matrix approximation for the two- 
particle Green's functions. This approximation scheme incorporates the resummation of the 
diagrams necessary in the presence of a strong repulsive short range component in nuclear 



forces and at the same time is thermodynamically consistent 



3,Q. 



In the zero tempera- 



ture limit, the in-medium T-matrix method gives results similar to the variational and BHF 

approaches 0, 0, Q, 13 

for the binding energy. At finite temperatures, however, it has 
the important advantage of automatically fulfilling the thermodynamic consistency relations, 
guaranteeing an unambiguous estimate of the thermodynamic quantities 15] . It has been 



employed in several works including studies of the modification o 



medium 



20, 



2l| . effective scattering 22|, pairing correlations |23l. 



nucleon properties in the 



sition 25J and nuclear matter binding and thermodynamics 




241, liquid -gas phase tran 



28|,129|, 
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Existing self-consistent T-matrix calculations take into account only the two-body in- 
teractions, with the exception of a first attempt using a purely phenomenological Lagaris- 
Pandharipande TBF 31|. In the present work we implement the self consistent T-matrix 
scheme with two different realistic NN potentials (CD-Bonn and Nijmegen) plus Urbana 
three-nucleon forces, which consist in an attractive part, based on the two-pion exchange, 
and a repulsive term related to the presence of a Roper resonance as an intermediate state. 

Within this approach we calculate the binding energy per particle for both symmetric 
and pure neutron matter, and give an estimate of the nuclear symmetry energy. Moreover 
we study the nucleon spectral function, self-energies and effective mass in comparison with 
the case of two-body forces only. 

The paper is organized as follows. The next section is a brief review of the Green's 
function formalism and the T-matrix approximation. In the third section we illustrate how 
the three-body forces are included in the T-matrix scheme, and in section IV we present 
our results on the nuclear EOS and the single particle properties, placing emphasis on the 
modifications due to the introduction of TBF. In the last section we end with a discussion 
and a short summary. 

At the moment we have considered the zero temperature case only. We plan to perform 
calculations at finite temperature, yielding thermodynamic quantities and a more complete 
nuclear EOS including TBF. This will be the object of a future publication. 



II. SELF-CONSISTENT T-MATRIX APPROACH 



In the diagrammatic expansion of the nucleon self-energy, the summation of the ladder 



dragrams at all orders leads to the T-matnx or ladder approximation Q 123 . The imaginary 
part of the self-energy is obtained as follows (for more details see |30[ and references therein) 

d 3 k dJ 
(2tt) 3 2tt 

-Im ((p - k)/2|T*(p + k, u + o/)|(k - p)/2)l \b(u> + J) + f{J)] A(k, J) 



ImE^p,^) 



Im <(p - k)/2|T K (p + k, u + o/)|(p - k)/2> 



Here b(u) and f(uj) are respectively the Bose-Einstein and the Fermi-Dirac distributions, 
A(p, uj) is the nucleon spectral function related to the single particle propagator through 



g<(p, w ) = /H4p )W ). 



(2) 
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The self-consistent in-medium two-particle scattering matrix T is defined as 

(k\T R W(P,Q)\\s!) = V(k,k') 
+ / ^^^(k,p)(p|Gr i?(A) (P,^)|q)(q|T i?(A) (P^)|k') , (3) 

where V(k, k') is the chosen nucleon-nucleon potential and the uncorrelated two-particle 
Green's function GVf is the product of two dressed one-body propagators 

(k'|G2 C<(>) (P,fi)|k) = 
i(2nf 5(k - k') / ^G <(>) (P/2 + k,Q- J) G <(>) (P/2 - k, u') . (4) 

J Z7T 

The single particle propagators are dressed via the Dyson equation 

G^)- 1 ( P ^) = c-|^-S^)(p^). (5) 

2m 

Formulae ([T]), and §5§ represent a closed set of equations which has to be solved with an 
iterative algorithm, until self-consistency is achieved. During the iteration the Fermi energy 
/i is adjusted to get the desired density p from 

Once stable solutions are obtained, the energy of the many-body system is calculated directly 



from the expectation value of the Hamiltonian (see |30|), avoiding the use of the Galitskii- 
Koltun sum rule. This is an important point because the sum rule, besides showing some 
dependence on the chosen energy integration cutoff, loses its validity in the presence of 
three-body forces. 



III. THREE-BODY FORCES 



The profound reason for the appearance of three-nucleon forces is that protons and neu- 
trons are not unstructured particles, as they are considered when constructing the various 
NN potentials in meson-nucleon field theory. When a nucleon interacts its internal structure 
indeed can be modified; this implies that there are interactions which are not anymore sim- 
ply of the nucleon-nucleon type. This is usually described via the excitation of one nucleon 
to another state or resonance, giving rise to processes which cannot be decomposed into 
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a sequence of NN diagrams. At low densities their contribution is negligible, while they 



34]. 



become relevant at densities involved in nuclear matter calculations 

All calculations aiming at realistic predictions of nuclear matter properties around and 
above the saturation point include TBF. There exist few works which implement TBF in 
nuclear calculations and test their effects on binding energy, pairing and single-particle 



pro per ties in different approaches, both in finite systems and in nuclear matter [11, 12, 13, 



35 



36, 



37 



38 



39, 



40] 



We consider here the approach of the Urbana three-nucleon potential developed by Carl- 



son et al 



3|. It is composed of two terms 



Vi 



(7) 



The first part is constructed from two-pion exchange, where a A resonance appears as inter- 
mediate state. Its contribution is attractive and dominates at low densities (below nuclear 
saturation). The second term, whose nature is more phenomenological, provides the repul- 
sion needed at higher densities. It is introduced to represent the contribution of several other 
diagrams, possibly including in an effective way also four- and more nucleon interactions. 



Alternative implementations and microscopic derivations of TBF are also available 
In more details, the Urbana TBF have the following structure 



34| 



2tt 
ijk 



A ( i X V' X 3k} {n ■ Tj, Tj-T k } + - [Xij, X jk ] [n- Tj,Tj • T k ] 

eye ^ ^ 



where 



and 



Yin 



(Tj + '/'(/•,, ) S tJ : 



^ijk 



jk) 



(8) 



(9) 



(10) 



eye 



The radial functions Y(r) and T(r) are the Yukawa and tensor functions respectively, the 
tensor operator is defined as Sij = 3 (<7j • ry)(<Tj • r^) — crj • crj, where is the unit vector 
of the distance between particles % and j. To determine the overall strength of the TBF 
and the relative strength between the two terms two parameters are present (A < and 
U > 0), to be tuned to reproduce the saturation properties of symmetric nuclear matter. 
Since different NN potentials lead to different saturation curves one should expect these 
parameters to depend on the particular choice of the two-body force. 
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The three-body interaction depends on the spatial, spin and isospin coordinates of the 
three nucleons, and in such a form cannot be used in the calculations. We then need to 
introduce some approximation and derive an effective two-particle potential. This can be 
done by averaging the action of the third nucleon, resulting in a mean field felt by the other 
two: 

K 3 //(<Z, <?') = E / < k ) V ^ 9. «0 , (11) 
where V 3 (k, q, q') is the Fourier transformed form of (JTj) and 

n(k)= f^G<(k,uj) (12) 

J 2,71 

is the particle momentum distribution. The sum over spin and isospin degrees of freedom 
just reminds us that V 3 has a non trivial structure in the a and r spaces which has to be 
taken care of (we didn't write explicitly spin and isospin indices). 

This average has to be performed for each of the three nucleons and over all their possible 
permutations, resulting in nine different terms. One has to pay particular attention to the 
spin-isospin and tensor dependence of the various averages and finally get, for each of the 
nine permutations, an effective potential of the form 

V?,/(9, <?') = V s R (q, q') + V?(q, q') (13) 
+ V 2 ?(q, q')* <T'T T' + Vg{q, q') S(q, q') r ■ r' , 

where V S R , V 2n , V 2 ? and Vg* are now scalar functions. 

Once we have obtained V^f (density dependent) we add it to the two-body potential in 

© 

V^V' = V + V e 3 ff . (14) 
and perform the T-matrix iteration. 



IV. BINDING ENERGY AND SINGLE PARTICLE PROPERTIES 

We perform calculations with two different parameterizations of the NN interaction, the 
CD-Bonn and the Nijmegen potentials. For both of them we compute the energy per particle 
directly from the expectation value of the interaction Hamiltonian, for symmetric and for 
pure neutron matter, with and without TBF. In the case of three-body forces we have tuned 
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the parameters A and U in (181) and (TTOT) in the symmetric case in order to reproduce the 
saturation density po and binding energy Eq. Since the averaging over the third nucleon in 
TBF terms represents a rather crude approximation, the resulting numerical values of the 
parameters of the TBF are different than in other approaches. 



A. Symmetric nuclear matter 



The energy per particle as a function of density for symmetric nuclear matter is shown in 
Fig. CD The calculations with only two-body forces fail to reproduce the correct saturation 
behavior, predicting a saturation density p = 1.47 po in the case of the Nijmegen potential 
and p = 1.79 po for CD-Bonn. After the inclusion of three- nucleon interactions the situation 
is significantly improved, with both curves saturatin g ar ound the phenomenological value 
Po = 0.16 fm~ 3 and yielding a correct binding energy 52j (Nijmegen Eg = —16.4 MeV and 
CD-Bonn Eb = —16.3 MeV). We notice that the effects of TBF are almost negligible at low 
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FIG. 1: (Color online) Energy per particle in symmetric nuclear matter as a function of density 



in units of the nuclear saturation density po = 0.16 fm ). T-matrix calculations are compared 



to the variational and BHF [l3( approaches, both including TBF. 



n 



densities and become more and more significant as the density increases. 

We plot for comparison also the results obtained with the variational method by Akmal 



et al. [5( and with the BHF approach by Baldo and Maieron 



131 ]. both using the Argonne 



vl8 potential. The overall agreement between all the calculations is good. 
When we compute, however, the incompressibility at saturation 

d 2 E\ 



K = 9\p 



dp 2 



(15) 



P=PO 



we get for both potentials Kq ~ 150 MeV, a value smaller than the usual estimates Kq = 
210 ± 30 MeV extracted from the experimental data and alternative calculation schemes. 

The structure (TSJ) of the one-particle propagator allows us to investigate the correlations 
induced in the dense medium by NN and three-body forces. These correlations are reflected 
in the energy dependence of the spectral function A(p,u), which in the presence of the 
short-range nuclear interactions differs substantially from the ^-function behavior that char- 
acterizes the quasiparticle approximation. In Fig. [2] we show the spectral functions when 
only two-body forces are present and when TBF are added for the CD-Bonn potential at 
three different densities and zero momentum. 
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FIG. 2: (Color online) Spectral function at zero momentum for CD-Bonn interaction and symmetric 
nuclear matter, at po, 2 po and 3pQ. 



The inclusion of TBF results in a further broadening of the peak of the spectral function, 
with a larger effect at high densities. This behavior reflects an increase of the scattering 
between the nucleons and it is common to all densities and all momenta. In Fig. [3] two 
cases with non-zero momenta are presented at the density 3 po- in the left panel (jp close to 
the Fermi momentum pp), the strongest reduction of the peak shows up. The modifications 
caused by TBF are then slowly disappearing for higher momenta, as illustrated in the right 
panel. The effect of the TBF on the spectral function is minor for p < p . Consequently, TBF 
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FIG. 3: (Color online) Spectral function at p = 340 MeV (left panel) and p = 510 MeV (right 
panel) for CD-Bonn interaction in symmetric nuclear matter. 



and their isospin dependence cannot explain the 
extracted proton spectral functions 
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strong dependence of the experimentally 
45J on the target nucleus. 
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FIG. 4: (Color online) Nucleon self-energy for CD-Bonn interaction in symmetric nuclear matter 
at po (left panel) and 3 po (right panel) 



We have analyzed the real part of the self-energy, which determines the shift with respect 
to the free dispersion relation 

p 2 

u p = — + Re E(p, Up) . (16) 
It is the sum of the Hartree-Fock self-energy and a dispersive contribution obtained via the 
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imaginary part (pQ) 

ReT,(p,uj) = J: HF (p) +V — . (17) 

J 7T LO — LO' 

The introduction of TBF does not have a big impact on the self-energy. As shown in Fig. H] 
(left panel) for low densities the modification is smaller than 5%, and it grows up to 10% at 
p = 3po (right panel). At all densities the effect is relevant for low momenta and basically 
vanishes above the Fermi surface. As expected TBF have a repulsive contribution, shifting 
the dispersive part of Re£(p, uj) up to less negative values. The behavior is completely 
analogous if the Nijmegen potential is employed. 
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FIG. 5: (Color online) Effective mass (in units of the nucleon mass m = 939 MeV) for the CD-Bonn 
interaction and symmetric nuclear matter at po (left panel) and 3 po (right panel) . 
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Finally, the effective mass is evaluated in the presence of TBF. In general it is obtained 
at each momentum from 

dp 2 2m* ^ ^ 

The effective mass as function of the momentum is displayed in Fig. [5] for the CD-Bonn 
potential at two different densities. As for the self-energy, the effect of the TBF is not 
dramatic. At po TBF smooth out the peak formed around Fermi momentum pp = 263 
MeV, but do not change the value at the Fermi surface m* m 0.85 m which is in agreement 



with the typical estimates from experiments 46J and BHF calculation with rearrangement 



terms included. Both with and without TBF the effective mass at 3 po gets substantially 
lower, with m* w 0.6 m at the Fermi momentum. Similar behaviors are observed for the 
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Nijmegen NN interaction with and without three-body forces: at po we found m* ~ 0.9 m 
and at 3 po we estimate m* ps 0.65 m . 



B. Neutron matter 



In Fig. [6] we show the results for the energy per particle in the case of a system composed 
of neutrons only. For the two choices of the NN potential, curves without and with TBF are 
displayed. When TBF are used, the parameters A and U are fixed by the calculation in the 
isospin symmetric case. The introduction of TBF does not change qualitatively the density 
dependence of E/A but makes the EOS more stiff. As in the symmetric case the results are 
compared to variational and BHF calculations (both including TBF). We find agreement 
at low densities, with the T-matrix results getting stiffer at very high densities. A similar 
effect, a stiffer result from the T-matrix than from the BHF or variational calculations in 
neutron matter has been noticed when using two-body forces only 17]. 
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FIG. 6: (Color online) Energy per particle in neutron matter as a function of density (in units of 
the nuclear saturation density po = 0.16 fm~ 3 ) for different potentials. The T-matrix calculations 



are compared to the variational [5] and BHF 



IS] 



131 ] approaches, both including TBF. 



An important quantity in determining the equation of state of isospin asymmetric nuclear 
matter is the symmetry energy. In particular its behavior is crucial for neutron star physics, 
since it affects the stability of neutron star matter and controls the position of the crust-core 
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transition 47|. The symmetry energy is defined from the energy per nucleon E/A as follows 



, . 1 (d 2 ElA\ . . 

where 5 = (p n — p p )/p is the so called asymmetry parameter. We calculate E sym via the 
parabolic approximation, which is usually employed to estimate the energy per particle for 



arbitrary asymmetries (whose reliability has been checked by Bombaci and Lombardo 481 ]) 



j(n,5) = j(n,5 = 0) + 5 2 E sym (n) . (20) 

The experimental measurement of the symmetry energy at saturation E sym (po) = So is not 
simple and still has a rather large uncertainty. A common estimate is So = 32 ±6 MeV 49]. 
The density dependence of the symmetry energy is of crucial importance for neutron star 
physics. It is fairly unknown, with large differences, specially at high densities, between the 



various models 



50] 



In Fig. [7] our result for the symmetry energy is shown in comparison with the variational 
and BHF approaches. When employing the CD-Bonn potential we find So = 39.7 MeV, for 
Nijmegen So = 37.1 MeV, both slightly higher than other estimates. With two-body forces 
only we find Sg ~ 30 — 32 MeV compared to 29 MeV from an earlier self-consistent T-matrix 



calculation 511 ] . 



Single-particle properties in neutron matter are modified by the TBF. We present the 
results obtained with the CD-Bonn potential, premising that the same considerations are 
valid for the Nijmegen NN interaction. In Fig. [8] the spectral function at zero momentum is 
shown for three different densities. We see that in pure neutron matter the introduction of 
three-body forces causes an enhancement of the quasi-particle peak of the spectral function 
at all densities p > po- The suppression of correlations is more evident at low densities and 
decreases as we go up to 3 p . 

For all momenta up to pp there is a narrowing of the peak. The enhancement of the 
quasiparticle strength can be quantified as 20 — 30%. In Fig. [9] the spectral functions for 
two other momenta are displayed with qualitatively similar effects. 

The suppression of the effects of TBF above the Fermi surface is evident when we consider 
the self-energy (1171) . plotted as a function of momentum in Fig. [H]] for p = po (which 
corresponds to pp = 331 MeV) and p = 3 p {j>f = 477 MeV). The self-energy is shifted 
down by 5 — 10 MeV by the TBF, i.e. the single-particle potential becomes slightly more 

13 
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Color online) Symmetry energy as a function of the density. The T-matrix results are 



compared to the variational [5[] and BHF [13( approaches. 
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FIG. 8: (Color online) Spectral function in neutron matter at zero momentum for the CD-Bonn 
interaction with and without TBF. 



attractive. This causes a reduction of the neutron effective mass (see Fig. [TTI) for momenta 
below pf- Similarly to the symmetric case, for po after the introduction of TBF the structure 
around the Fermi momentum disappears, and the value at the Fermi surface is shifted from 
m* ~ 0.9 to m* « 0.85 . At 3 po there are only small differences when including TBF. With 
or without TBF the effective mass at the Fermi surface is found to be m* w 0.75. 
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FIG. 9: (Color online) Nucleon spectral function in neutron matter at 3po for p = 280 MeV (left 
panel) and p = 340 MeV (right panel). 
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FIG. 10: (Color online) Different contributions to the self-energy ReS(p, uS) at po and 3po in 
neutron matter (CD-Bonn interaction). 



V. SUMMARY AND CONCLUSIONS 



The equation of state of symmetric and asymmetric nuclear matter is of crucial interest 
in the study of heavy ion collisions and the modeling of neutron stars. The self-consistent 
Green's function technique, in which the strong correlations induced by the nuclear interac- 
tions are resummed in the in-medium T-matrix, accounts for the modified particle properties 
in the dense medium and provides consistently the macroscopic thermodynamic observables. 
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FIG. 11: (Color online) Effective mass in neutron matter (in units of the nucleon mass m = 939 
MeV) for the CD-Bonn interaction at po (left panel) and 3 po (right panel) . 



Within this approach, we calculate the properties of isospin symmetric and asymmet- 
ric nuclear matter at zero temperature, focusing on the impact of three-body forces on 
the density dependence of the energy per particle and of the single-particle properties. 
We employ two different NN potentials, CD-Bonn and Nijmegen, and the Urbana semi- 
phenomenological TBF. The results are qualitatively independent of the starting NN poten- 
tial, so that the present conclusions are valid for both of them. In table [I] we summarize the 





Po 
[fm^ 3 ] 


[MeV] 


K 
[MeV] 


So 
[MeV] 


m* /m 


experiment 


0.16 ±0.01 


-16 ±1 


210 ± 30 


32 ±6 


» 0.8 


CD-Bonn 


0.287 


-19.9 


70 


32.2 


0.90 


CD-Bonn + TBF 


0.171 


-16.3 


148 


39.7 


0.87 


Nijmegen 


0.235 


-18.4 


76 


30.5 


0.87 


Nijmegen + TBF 


0.164 


-16.4 


158 


37.1 


0.90 



TABLE I: Summary of the main features of the nuclear matter EOS obtained with and without 
TBF, compared with experimental estimates taken from 49]. 



most important parameters of the equation of state calculated with and without three-body 
forces. There is an evident improvement after the introduction of TBF, in particular in the 



1(3 



description of the saturation properties of symmetric nuclear matter. The incompressibility 
at saturation Kq is perhaps the only unsatisfactory result, being smaller than the current 
estimates. We have tested also a stiffer set of TBF parameters but we are not able to in- 
crease K Q without spoiling the neutron matter equation of state, which has already a stiff 
slope. 

In general, in a dense medium three-nucleon forces play a significant role. The modifica- 
tions of the EOS and of the single-particle properties are larger when the density increases. 
Analyzing the momentum dependence we notice that the effects of TBF are maximal around 
the Fermi surface, and are strongly suppressed for high momenta. 

At high densities the nucleon spectral function is modified significantly. This can change 
the estimate of the in-medium nucleon-nucleon cross section, which is relevant for the de- 
scription of heavy-ion reactions and the cooling of neutron stars. The effects depend on 
the isospin asymmetry of the system. It is a broadening of the quasi-particle peak of the 
spectral function in the case of symmetric matter, which reflects the increase of scattering 
between particles. On the contrary for neutron matter the peak is enhanced suggesting a 
stronger quasiparticle behavior. At subsaturation densities, relevant for finite nuclei, the 
effect of TBF on the spectral function is small. Thus, TBF cannot explain the observed 



strong dependence of t 
neutron ratio in nuclei 



re proton spectral functions on the average density and proton to 



4l|. 



We also analyze the dispersive self-energy and the resulting optical potential is more 
repulsive (attractive) with TBF in symmetric (neutron) matter. The effective mass is not 
very sensitive to the effects of TBF. Its momentum dependence is slightly smoothed out 
with the introduction of TBF. 

All the present results are at zero temperature. Within the T-matrix approach, it is 
possible to investigate the thermodynamic properties of such systems at finite temperature, 
including TBF forces into the calculations. TBF are a necessary ingredient in order to reach 
a realistic EOS of dense and hot nuclear matter. 
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